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We express the optimization of entanglement witnesses for arbitrary bipartite states in terms 
of a class of convex optimization problems known as Robust Semidefinite Programs (RSDP). We 
propose, using well known properties of RSDP, several new sufficient tests for separability of mixed 
states. Our results are then generalized to multipartite density operators. 
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I. INTRODUCTION 

Entanglement, first noticed by Einstein, Podolsky, 
and Rosen Q, is at the heart of quantum mechanics. 
Quantum teleportation, superdense coding and cryptog- 
raphy J} are achieved only when one deals with insepa- 
rable states. Thus, the determination and quantification 
of entanglement in a composite quantum state is one of 
the most important tasks of quantum information the- 
ory. In the past years a great deal of effort have been 
made in order to obtain the characterization of separable 
bipartite mixed states Q . A finite-dimensional bipartite 
density operator pab € B(Ha <S> Hb) (the Hilbert space 
of bounded operators acting on Ha ® Hb) is separable 
iff it can be written as a convex sum of separable pure 
states: 

Pab = y]pi\ipi) AAjtpi] ® \<t>i) BB{<t>i\ (1) 

where {pi} is a probability distribution and \ipi}A, \<Pi)B 
are vectors belonging to Hilbert spaces Ha and Hb, re- 
spectively. Despite the simplicity of this definition, none 
operational necessary and sufficient criterion have been 
found for the separability problem until now. More- 
over, it was showed by Gurvits that this problem 
is NP-HARD. Therefore, we should not expect to find 
a polynomial-time algorithm which determines for any 
state pab, with certainty, if it is possible to decompose 
it in the form of equation (1). 

A particularly useful concept is that of entanglement 
witness (EW). According to an operator pab is 
entangled iff there exists a self-adjoint operator W S 
B(Ha ® Hb) which detects its entanglement, i.e., such 
that TriyVpAB) < and Tr(Wa AB ) > for all <t A b sep- 
arable. This condition follows from the fact that the set 
of separable states is convex and closed in B(Ha <8 H B )- 
Therefore, as a conclusion of the Hahn-Banach theorem, 
for all entangled states there is a linear functional which 
separates it from this set. Unfortunately, it is not known 
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how to construct EW in a canonical way and in polyno- 
mial time for every entangled state. Actually, since such 
method would solve the separability problem, it can not 
exist as long as the strong conjecture P ^ NP is true. 

In this paper, we show that the search of EW for ar- 
bitrary mixed states is indeed NP-HARD. We introduce, 
in the context of quantum information, a class of con- 
vex optimization problems known as Robust Semidefi- 
nite Programs (RSDP), whose NP-hardness in most of 
the cases was already proved. This family generalizes 
the important Semidefinite Programs (SDP), which have 
been increasingly used in quantum information problems 
(see, for example, and |To|')- 

The paper is organized as follows. In section 2 we 
briefly recall the definition of EW and define our concept 
of optimal entanglement witness; in section 3 we state 
the basic facts about Robust Semidefinite Programs, ex- 
press the optimization of EW as a RSDP and provide 
a first approximation in terms of SDP for the problem, 
which yields a new sufficient criterion of separability; in 
section 4 we parametrize all possible approximations of 
the RSDP in terms of a multiplier matrix, reducing the 
search space of approximation scenarios; in section 5 we 
generalize our results to multipartite states. Finally, in 
section 6, we present our conclusions and suggest direc- 
tions for further research. 



II. OPTIMAL ENTANGLEMENT WITNESS 

A hermitian operator W £ B{Ha ® Hb) is an entan- 
glement witness if |5j [6|: 

1. a{iP\ ®b {<j)\W\4>) B ® W)A > 0, for al l states \ip) S 
H A and \<f>) e H B . 

2. W has at least one negative eigenvalue. 

3. Tr(W) = 1. 

Condition one assures that Tr(WaAB) > for all separa- 
ble states cj ab • Condition two implies that Tr(WP) < 
at least for one entangled state, for example, the projec- 
tor on the eigenspace associated with the negative eigen- 
value. The third condition is important in order to com- 
pare different EW. 
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Definition 1 A hermitian operator W PAB is an optimal 
EW for the density operator pab if 

Tr(W PAB p AB )<Tr(Wp AB ) (2) 

for every EW W . 

Although the above definition of OEW is different from 
the one introduced in pfl ,the optimal EW of both criteria 
are equal. According to 0, W is optimal iff for all P > 0, 
W = (1 + e)W - eP is not an EW. 

III. ROBUST SEMIDEFINITE PROGRAMS 

In this section we will express the search of an opti- 
mal EW for an arbitrary state pab in terms of a robust 
scmidefinite program (RSDP). A semidefinite program 
(SDP) consists of minimizing a linear objective under a 
linear matrix inequality (LMI) constraint, precisely, 

minimize c^x subject to 

m 

F(x) =F„+J2 x iFi > ( 3 ) 

i=l 

where c € C m and the hermitian matrices Fj = Fj 6 c nxn 
are given and x € C m is the vector of optimization vari- 
ables. F(x) > means F(x) is hermitian and positive 
semidefinite. SDPs are global convex optimization pro- 
grams and can be solved in polynomial time with interior- 
point algorithms For instance, if there are m opti- 
mization variables and F(x) is a nxn matrix, the number 
of operations scales with problem size as 0(m 2 n 2 ). SDPs 
have already been used in different problems of quantum 
information theory [t| and also in the separability prob- 
lem |l0l| . An important generalization of (3) is when the 
data matrices Fi are not constant, i.e., they depend of a 
parameter which varies within a certain subspace. This 
family of problems, known as robust semidefinite pro- 
grams, is given by: 

minimize c^x subject to 

m 

F(x, A) = F (A) + XiFi(A) > 0, VA e V (4) 

i=l 

where I? is a given vectorial (sub)space. Note that prob- 
lem (4) is more difficult to solve than (3) , since one must 
find an optimization vector x such that F(x, A) is pos- 
itive semidefinite for all A e T>. One often encounters 
SDP in which the variables are matrices and in which 
the inequality depends affincly on those matrices. These 
problems can be readily put into the form (3) by intro- 
ducing a base of hermitian matrices to each matrix vari- 
able. However, since most of optimization solvers [TT| ad- 
mit declaration of problems in this most general form, it 
is not necessary to write out the LMI explicitly as (3) , but 



instead make clear which matrices are variables. Equal- 
ities constraints involving the optimization variables can 
also appear in (3) and (4) without any further computa- 
tional effort. We can now enunciate the main result of 
this paper. 

Theorem 1 A state pab € B(Ha <8> Hb) is entangled, 
i.e., can not be decomposed as (1) iff the optimal value of 
the following RSDP is negative: 

minimize Tr(WpAB) subject to 

d.A d A 

Y^L a i a i W H ^ °> Tr(W) = 1, for all a t e C (5) 

»=i 3=1 

where dA is the dimension of Ha, W%j —a (i\W\j)A € 
B{Hb) and \ j) a is an orthonormal base of Ha- If Pab 
is entangled, the matrix W which minimizes Tr(WpAB) 
is the OEW for pab- 

proof: First we have to show that (5) is a genuine RSDP. 
Note that Wij =a (i\W\j)A and the objective Tr(W pab) 
are both linear in the matrix variable W. Thus, (5) 
can be put into the form (4), where T> in this case is 
C dA . We know that a state pab is entangled iff there 
exists an operator W such that Tr(WpAB) < and 
a{*1>\ ®b {<f\W\<j>) B ® \^)a > for all states \ip) e H A 
and \<j>) £ H B - Therefore, the matrix a(V'|W / I'0)a has 
to be semidefinite positive for all \i/j)a G Ha- Letting 
\iP)a — J2j a*j\j)A, where \j)a is an orthonormal base of 
Ha, it is straightforward to show that the optimal W 
given by (5) is the OEW of pab ■ QED. 

In spite of the similarity between (3) and (4), RSDPs 
are in general very hard optimization problems. Actually, 
it was proved that robust semidefinite programs like (5) 
are NP-HARD [l^. 

Corollary 1 The determination of the OEW for an ar- 
bitrary state pab is a NP-HARD problem. 

Since (5) is computationally intractable, it is natural to 
search for approximations of it in terms of SDPs, which 
are very efficiently solved. These relaxations of RSDP 
have been intensively studied ^| m the past years and 
can be classified as probabilistic or deterministic. In 
this paper we will focus in the deterministic relaxations, 
where (4) is replaced by an inner convex approximation 
described by a linear matrix inequality constraint. This 
inner approximation is then used to find an upper bound 
to the optimal value of (4). The probabilistic approach, 
which yields outstanding results on the separability prob- 
lem, will be reported elsewhere. As a first example of 
such relaxations, consider the following adaptation of 

0: 

Theorem 2 A density operator pab is entangled and the 
optimal value of W is an EW for it if the result of the 
following SDP is negative: 

minimize Tr(W pab) subject to (6) 
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1- W kk >0, k = l,2,..,d A 

2 - di=i Wkk ± *r(- w kj + w 3k) > 0, l < k + j < d A . 

3- dJ=T w kk ± Tf(W kj - W jk ) > 0, l<k^j<d A . 
proof: Note that: 

1 



^2a* k a j W k j= 

k,j i<k<j<d A 



+ a%ajW kj + a*a k W jk + 



d A -l 
1 



\a k rW k k 



d A -l 

Thus, a sufficient condition to (5) is: 

l&le) (e\W kj \e) 



\a k \ z W kk 



> 



M{\e)) 



{e\W jk \e 



> 0, V|e) G H B 



This matrix M(\e}) is positive semidefinite iff its diagonal 
entries and determinant are greater than or equal to zero. 
From condition 1 it follows that Mu > and M22 > 0. 
From conditions 2 and 3: 

(e| rf7^T |e>(e| rfI^T |e> " 2 l max {{e\% W kj+W jk )\e), 

(e\^§(W kj - W jk )\e}}} 2 > (e\W kj \e)(e\W jk \e) 
det(M(\e)) > 0. QED. 

We present now the first example of our methodology. 
We used MATLAB and the package SEDUM1 fllj to im- 
plement and solve the SDP. 



A. Bell State 

We consider the following Bell state |\&) = (|0) ® |0) + 
|1) <g> |l))/\/2. It is well known that (#|Wj*>(*||#) = 
-1/2 where W\ 9)W = — 1(|00)<11| + |11)<00|) + 
i(|0f)(0f| + |10> <10|) is the OEW for this state. Solving 
the SDP of theorem 2, the following EW was found 



W = 



0.1057 -0.2887 

0.3943 

0.3943 

-0.2887 0.1057 



Since program (6) is only a relaxation of (5), W is not 
the OEW for |*>, (*|W OJ)t |#) = -0.1835. 



OAS 



1/3 





2/3 



The positive partial transpose criterion show that p A s is 
entangled, while a A B is separable. Using the SDP (6), 
we have found the following EW for p A B- 



W 



0.1752 

0.1752 -0.2478 

-0.2478 0.0513 

0.5982 



(7) 



where Tr(Wp AB ) = -0.0313. The method has also 
succeeded in the state <j A b, as the optimal value for 
Tr(Wa AB ) founded was 2.7330a;10- 5 . 



IV. COMPLETE FAMILY OF PARAMETRIZED 
RELAXATIONS 

It must be stressed that theorem 2 is only one of the 
possible approximations of (5). In fact, every relaxation 
of the RSDP constitutes a different method of EW con- 
struction and, therefore, a new sufficient criterion of sep- 
arability. In this section we will show that all these possi- 
ble relaxations can be parametrized in terms of a family 
of matrices. However, in order to provide such method, 
it is necessary firstly to introduce some standard results 
concerning robust semidefinite programs. One partic- 
ularly important representation of robust linear matrix 
inequalities is the Linear Fraction Representation (LFR) 
[l3| |T^ . It was showed that every matrix F £ C nxc 
which depends rationally of a varying parameter S G C k 
can be expressed as [T^ |: 



F{5) = A + BA(I - DA) _1 C 



(8) 



where A G C nxc , B G C nxN ,C G C Nxc and D G C NxN 
are constant matrices, n, r k and N = n + ... + r k are 
integer numbers and A is the following diagonal matrix: 

A = Diag(5iI ri ,...,5 k I rk ) 

We can now express problem (5) in terms of a LFR. 

Theorem 3 A state p A s is entangled and the optimal 
value ofW is the OEW for it iff the result of the following 
RSDP is negative: 

minimize Tr(Wp A B) subject to 



B. Isospectral States 



We now consider the two isospectral matrices: 



p A B 



1/3 

1/3 1/3 

1/3 1/3 





F(A) = BA(I — DA)~ X C > 0, VA G V (9) 

where 

B = [W n - W nl W 12 ... W nn dA , dAdB ] (10) 



C 



dA ,d%d A L®I dA 



it 



(11) 
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D 



^d A d\-d A d B L^®I dA d E 

®d A d B ;d A d B Od A d B ;d A d E 



(12) 



a *d A Id A) 



A = Diag (a 1 I dAdB , a dA Id A d B , a\I dA1 

' (13) 

and T> is the subspace of diagonal matrices in the form 
of (12), where dj € C. O p - q and I p stand for thepxq zero 
matrix and the pxp identity matrix, respectively. L is an 
auxiliary matrix given by: 

L = [i,i,...,i] e h b 

proof: We must show that the LMI of (8) is equivalent 
to the LMI of (5). In order to do that, we will use con- 
structive formulas of addition and multiplication of LFRs 
presented in the appendix. Each quadratic term from the 
LMI of (5) can be written as 

a*ajWij = {0 + 1 x a* (I - x a*) -1 x l}x 



{0 + Wij xa*{I-0x a*)- 1 x /} 

A LFR to each term and then to the whole expression 
can now be obtained using the addition and the multi- 
plication formulas, respectively. QED. 

There are several approximations for robust linear ma- 
trix inequalities (RLMI) which are described as LFR 
[l3j . One of particular importance is the Full Block S- 
Procedure |17| : 

Theorem 4 (Full Block S-Procedure TF\ l) The matrix 
F(A) = A + BA(I - DA)~ l C is well posed and satisfies 



I 


t 


' X ' 




/ 


_F(A)_ 




X 




_F(A)_ 



< 0, VA e V (14) 



iff there exists a multiplier 
P = 



Q S 
R 



Theorem 5 A state pab is entangled iff there exists a 
multiplier matrix (15) such that (16) and (17) hold, with 
X = — I . The matrices appearing in the LMI (17) ore 
given by equations (10 — 12) and the matrix A is given 
by (13). 

proof: Noticing that F(5) > is equivalent to (14) if 
X = —I, the result follows easily from the application 
of the Full Block S-Procedure (theorem 4) in the RSDP 
(9). QED. 

The families of matrices P such that (16) is satisfied 
parametrize all possible relaxations of (5). Although the 
determination of all such matrices is not a trivial prob- 
lem, it is a lot easier than (14) and it is current subject 
of intensive research. Further information on possible 
choices of the matrix P can be found in [13. As an 
example, we consider now the most simple family of ma- 
trices P for which (16) holds. In quantum mechanics one 
usually deals with normalized states (ip\ip) = 1. There- 
fore, the matrix A satisfies A^A < /. In this case, the 
following matrix P gives an approximation of (5): 



P = 



-I 
/ 



(15) 



V. MULTIPARTITE ENTANGLEMENT 

So far we have only considered the bipartite case. In 
this section we generalize the previous results to multi- 
partite states. A density operator pi...„ € B(Hi®...®H n ) 
is separable if it can be decomposed as: 

Pl..,n = y^Pi\lpi)ll(lpi\ ® •■■ <8> \lpi)nn(^i\ (18) 

1=1 

Since the set of multipartite mixed separable states is also 
convex, it is possible to apply the Hahn-Banach theorem 
and establish the concept of EW in a straightforward 
manner |18|. 

Theorem 6 A state p\... n E B(Hi (g> ... (g) H n ) is entan- 
gled, i.e., can not be decomposed as (18) iff the optimal 
value of the following RSDP is negative: 

minimize Tr(W pi... n ) subject to 



with 



A 
I 



i f 



P 



A 
I 



> o, va e v 



(16) 



EE- E E {<■ 



i„_i=l j„_i 



(19) 



such that 



> 



" / 


o ■ 


t r 


A 


B 







I 




. c 


D _ 





A 

A 

Q S 

R 





- 1 


" 




A 


B 







I 




. c 


D _ 



< 



(17) 



We can now express all possible relaxations of (5) in 
terms of the multiplier matrix P as follows: 



Tr(W) 



1, Va. ifc eC, l<k<n 



where d n is the dimension of H n , Wi 1 ,,,i rl _ 1 j 1 ...j n _ 1 = 
i(i|®...®„_i(i|W|j)„_i®...®|j)i € B(JTi(g)...®F„_i) 
and \j)k is an orthonormal base ofHk- If pi... n is entan- 
gled, the matrix W which minimizes Tr(W p\ „) is the 
OEW for Pl ... n . 
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proof: We know that a state pab is entangled iff there 
exists an operator W such that Tr(Wpi... n ) < and 
1 (ip\®...(g) n (ip\W\il;) n (g)...(g)\ip) 1 > Ofor all states \ip) k S 
H K . Thus, the matrix i(V>| ® — ®n-i (V'I^IV')n-i ® •■■ ® 
(V')i > has to be semidefinite positive for all \ip)k G -Hif ■ 
Letting = 2j a jili)fc> where |j)fe is an orthonormal 
base of 7i & , it is straightforward to show that the optimal 
W given by (19) is the OEW of Pl ... n . QED. 

Relaxations for (19) can be obtained using the same ar- 
guments exposed before. Since the RLMI of (19) is poly- 
nomial in the varying parameters, it can be expressed as 
a LFR and we can apply the Full Block S-Procedure to 
the multipartite case. Therefore, all possible determinis- 
tic approximations of (19) can also be parametrized by 
the matrix P. Further results concerning the application 
of possible families of parameterizations of P in the op- 
timization of EW will be reported elsewhere. 



VII. APPENDIX: COMBINATION OF LFRS 

We provide in this appendix some simple combination 
rules for addition and multiplication of LFR used in this 
paper. Consider two matrix described by the LFR for- 
mat: 

Fi(S) = A + BiAi(I - D i A i ) -x Ci 
The sum of F± and F2 has the LFR: 

F(S) = A + BA(I — DA)~ 1 C (20) 

with 

A = Ax + A 2 , B = [B x B 2 ], C= [d C 2 ] t 



VI. CONCLUSION 



D = Diag(D 1 ,D 2 ), A = Diag(A u A 2 ) 



In this paper we have introduced, in the context of 
Quantum Information, a new class of optimization pro- 
grams (RSDP), showing that the determination of the 
OEW for an arbitrary state is NP-HARD. Several possi- 
ble deterministic approximation scenarios have been pro- 
posed to it, yielding new sufficient criteria of separability. 
Our results were then straightforwardly generalized to 
multipartite states. It was also showed that all sufficient 
criteria of separability might be parametrized by a matrix 
which satisfies a much simpler linear matrix inequality. 
Therefore, it is of great importance a systematic study of 
all possibles families of parameterizations for this matrix. 



The product of F\ and F 2 is given by (20) with 
A = A X A 2 , B = [B X A X B 2 ], C = [C X A 2 C 2 ] 



D = 



D Q l %f , A = Diag(A u A 2 ) 
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